What is the temperature dependence of the Casimir effect? 
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There has been recent criticism of our approach to the Casimir force between real metallic sur- 
faces at finite temperature, saying it is in conflict with the third law of thermodynamics and in 
5_l , contradiction with experiment. We show that these claims are unwarranted, and that our approach 

Oh has strong theoretical support, while the experimental situation is still unclear. 
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>* . I. INTRODUCTION 

IT) 

Recently there has been a comment l] that has strongly criticized our work Q on the temperature dependence of 
\Q \ the Casimir force between a spherical lens and a plate, both coated with real metals (Au for example). This follows 
a long paper by some of the same authors Q, which also criticizes our work as being in conflict with fundamental 
thermodynamical requirements, particularly the third law of thermodynamics or the Nernst heat theorem. Further- 
more, they claim that their recent experimental results are completely consistent with their theoretical approach, and 
in contrast completely rule out our theory. The purpose of this note is to respond to these criticisms, and emphasize 
Oh, that their theoretical claims are invalid, while the experimental situation is still too indecisive to draw definitive 
1 conclusions. 

First, let us make a statement about the physical nature of the controversy. The conventional approach, dating 
back to Lifshitz and reinforced by Schwinger et al. describes an ideal metal by a formal e — > oo limit, where e 
is the permittivity of the material. Mathematically, the limit is taken in such a way that both the TE and TM zero 
modes (that is, the m = term in the Matsubara sum) contribute. Recently, however, it has been recognized that 
this is not correct, and that the TE zero mode cannot contribute either for finite e or for a real metal. Consequently, 
in the ideal metal limit, this recognition gives rise to a linear temperature term in the pressure at low temperature, 
and a nonzero value of the entropy at zero temperature. However, real metals are not described by this ideal limit, and 
even though only one zero mode, the TM one, can be present, as one sees by either thermodynamic or electrodynamic 
considerations, the contradiction with the Nernst theorem is not present for actual conductors. 

In contrast, the authors of 0,13 propose two alternatives to describe the reflection coefficients that enter into the 
Lifshitz theory, the impedance approach and the plasma model. In the former the transverse momentum dependence 
in the surface impedance is simply disregarded (the so-called Leontovich approach), while in the second alternative 
the plasma dispersion relation with no relaxation is used for the dielectric constant. Neither of these omissions is in 
accordance with the properties of real materials, as we will detail in the following. 
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II. LIFSHITZ FORMULA 



We start by remarking that the Lifshitz formula for the force between two parallel dielectric slabs is essentially 
geometrical, beingdetermined entirely by multiple reflections at the interfaces between the parallel media (for example, 
see Sec. 3 of Ref. Q). Now there is no doubt that one can use surface impedances to describe the reflection coefficients 
appearing in the theory, for these are merely the linear relations between transverse electric and magnetic fields at 
the surface of the metal. It is well known that use of either bulk permittivity or surface impedance are equivalent 
in describing the reflectivity properties, provided that one incorporates transverse momentum, k±. The boundary 
conditions implied by Maxwell's equations provide the precise connection between these alternative descriptions 
In principle both bulk permittivity and surface impedance should depend on k±, although theory and optical data 
suggest that the dependence of permittivity on kj_ is rather small. This would then imply a substantial dependence 
on k± in the surface impedance. Indeed, in Ref. we show that the exact impedance also leads to zero contribution 
of the TE zero mode. 



A. Transverse momentum dependence of the surface impedance 

The authors of Ref. P, |j| ignore the transverse momentum dependence at zero frequency, because they believe 
that the "mass shell" condition (w/c) 2 — k\ > must be satisfied. However, this relation can only be valid for real 
photons. In computing the temperature dependence, one is evaluating Green's functions periodic in imaginary time, 
with period (3 = 1/kT |lfj [. The Matsubara frequencies, £ m = 2irm/ [3, being imaginary, evidently break the mass-shell 
condition. Physically, evanescent waves or virtual photons are responsible for the thermal Green's functions. The 
inclusion of spatial dispersion, as required by Maxwell's equations, has been considered in detail in recent publications 
[ill IT^ and the effects of the exclusion of the TE zero mode at both high and low temperatures are manifested. 



B. Plasma model 

The authors of Ref. 0, Q also state that the plasma model is equally good at describing the optical properties 
of real metals. The plasma model is a special case of the Drude model, with the relaxation parameter neglected. 
In fact, optical data clearly show that this limit is inadequate: The Drude model describes the permittivity of real 
metals up to £ « 2 x 10 15 rad/s, while the plasma model fails for frequencies below 5 x 10 13 rad/s. It is precisely 
the low-frequency part of the spectrum that is relevant for the discussion of the disputed zero-mode contribution. 
Those authors argue that the neglect of relaxation is appropriate because the zero temperature limit of the relaxation 
parameter for an ideal metal (the Bloch-Griineisen law) is zero. To this we reply that the relaxation parameter for 
real materials cannot vanish at zero temperature because of the scattering by impurities, and further that all present 
and nearly all proposed experiments are carried out at room temperature, so the temperature dependence of this 
parameter should be irrelevant. Moreover, it has recently been demonstrated [l"H that at sufficiently low temperature 
the residual value of the relaxation parameter does not play a role, as the frequency characteristic of the anomalous 
skin effect becomes dominant. In contrast to the ad hoc procedure advocated in Ref. 0,0, which in addition to 
using physically inadequate models, employs an unjustified extrapolation from the infrared region (that is, the plasma 
model, which disregards relaxation, is extrapolated from that region down to zero frequency), we use real optical 
data at room temperature for the bulk permittivity. Elsewhere it has been emphasized that a dielectric function valid 
over a wide frequency range and not merely the infrared must be employed, and indeed that frequencies very small 
compared to the characteristic frequency play a dominant role in the temperature dependence [Tl|. Of course, there 
is no data at zero frequency, but physical requirements (thermodynamics and Maxwell's equations) do not allow the 
TE zero mode to be included in the Lifshitz formula. 



III. THERMODYNAMICS 



As far as we can understand there is no disagreement about the validity of the Lifshitz formula for the Casimir 
force as long as the dielectric constant is finite. For this latter situation there is no m = or zero frequency TE 
(transverse electric) mode, in agreement with Maxwell's equations of electrodynamics, and inclusion of such a term 
would violate the third law of thermodynamics \H\ . 

The controversy arises in the limit of a perfect metal. It is well known that in this limit the mathematics of the 
limiting process can be ambiguous due to its singular nature. So the limit of a perfect metal can in one way be taken 
by letting e — > 00 when regarding the contribution from the m = term. Since this term for all finite e is zero, its 
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limiting value is also zero. However, the other way is to take the limit m — > while e = oo [5(- In the latter case a 
non-zero contribution arises. The opinion of Bezerra et al. 0] is that the latter procedure is the proper one while we 
have the opposite opinion. 

So the controversy left is for perfect metals of infinite extension. One of our arguments is that this special case is 
a limiting case of more realistic models with relaxation. As mentioned earlier, this latter limit can come in conflict 
with the third law of thermodynamics. But even in this limit the possible violation is ambiguous or not obvious. On 
one hand the entropy remains zero for temperature T = in the limit e — > oo. On the other hand taking the limit 
T — > with e — oo will result in entropy different from zero. The latter is connected to the increasing and diverging 
slope of the entropy function as e — > oo close to T = 0. 

However, for realistic models for metals with relaxation (or finite conductivity) this controversy is not present. As we 
have found earlier the entropy then goes smoothly towards zero as T — > 0. Bezerra et al. pj correctly point to the fact 
that the so-called M1M (modified ideal metal) violates the third law of thermodynamics. (This we stated explicitly 
in Ref. |l3j|.) This violation is then in the sense discussed above, i.e., it is somewhat ambiguous depending upon 
how the limiting process is performed. It is also correct that according to our approach real metals approach MIM 
for increasing separation as then lower frequencies become more important. (MIM is the limit where the dielectric 
constant e — > oo with the TE zero-mode equal to zero.) But the crucial difference between real metals and MIM is 
that the former includes relaxation by which there will be no violation of the third law of thermodynamics (at least 
for finite separation of metal plates and for infinite separation where there is no force or interaction energy left). 
Moreover, numerical studies have shown that the MIM model disagrees strongly at short distances with values for the 
Casimir force calculated for real metals even for zero temperature. 

Bezerra et al. [lj further claim that our unstated assumption is that perfect crystals with no defects or impurities do 
not or cannot exist, and, therefore, that thermodynamics can be violated for them. At best this is a misinterpretation 
of our arguments. There is no violation of the third law in this connection as quantization of lattice vibrations on a 
perfect lattice yields no problem; in fact it can be performed exactly. A reason for this is that the standard lattice has 
a finite dielectric constant and there is no ambiguity. However, if the perfect lattice turns into a perfect metal (with no 
relaxation) we can again sec the problem discussed above. Then there will be no well defined thermal equilibrium as 
any steady current can stay on forever. And we do not think that Nernst had this latter very special and unrealizable 
situation in mind when he formulated his theorem. 

Thus, we believe that the results obtained by excluding the TE zero mode are consistent with the third law of 
thermodynamics. The entropy at zero temperature vanishes, except in the limit of a perfect conductor of infinite 
extent, for which thermal equilibrium can never occur. There is a region in which the Casimir entropy is negative, 
but this is a rather familiar phenomenon, reflecting the fact that we are describing only part of the complete physical 
system. (Although Svetovoy and Esquivel an d Sernelius ^3 agree with our conclusion that the entropy vanishes 
at zero temperature for real metals with no TE zero mode, and thus there is no contradiction with the Nernst theorem, 
the former authors believe that there is a thermodynamic difficulty with the negative entropy region.) 



IV. THEORETICAL SUPPORT 

Two other recent papers also lend support to our point of view. Jancovici and Samaj 0] an d Buenzli and Martin 
|15| have examined the Casimir force between ideal-conductor walls with emphasis on the high-temperature limit. Not 
surprisingly, ideal inert boundary conditions are shown to be inadequate, and fluctuations within the walls, modeled 
by the classical Debye-Hiickel theory, determine the high temperature behavior. The linear in temperature behavior 
of the Casimir force is found to be reduced by a factor of two from the behavior predicted by an ideal metal. This 
is precisely the signal of the omission of the m = TE mode. Thus, it is very hard to see how the corresponding 
modification of the low-temperature behavior can be avoided. 

Further support for our conclusions can be found in the very recent paper of Sernelius 16], who calculates the van 
der Waals-Casimir force between gold plates using the Lindhard or random phase approximation dielectric function. 
Spatial dispersion plays a crucial role in his calculations. For large separation, the force is one-half that of the ideal 
metal, just as in the calculations in 0,0. For arbitrary separations between the plates, Sernelius' results numerically 
nearly exactly coincide with his earlier ones [6l Il7j where dissipation (i.e., nonzero relaxation parameter) is included 
but no spatial dispersion. In his new calculation, the inclusion of dissipation has negligible additional effect. His new 
results thus essentially coincide with ours. 
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V. EXPERIMENTAL CONSTRAINTS 

Thus, we claim that there is now overwhelming theoretical evidence that a TE zero- mode should not be included in 
the Lifshitz formula when describing the interaction between real metal surfaces. This makes the purported exclusion 
of this theory by recent experiments, as forcefully stated in Refs. 0,0], most difficult to understand. We believe that 
the resolution of this conundrum lies in an insufficient appreciation of the backgrounds making Casimir measurements 
so difficult. We and others have mentioned the difficulty in accurately determining the absolute sphere-plate 
separation. This may be especially so since the roughness of the surfaces is much larger than the precision stated in 
the determination of the separation. Nonlocality in the electric properties of the surfaces, as for example seen in the 
calculation of Rcf . means that the location of the surface cannot be specified to better than some effective shielding 
length. Interferometric methods of determining distance may also be somewhat uncertain because of penetration of the 
surface by light. Also accurate determination of a small difference between experimental values at room temperature 
and purely theoretical values at T = gives rise to further difficulties. In short, fitting data precisely to a preferred 
theory wherein various unknown experimental parameters must be determined does not constitute decisive evidence 
in favor of that theory; nor does the poorness of fit to an alternative theory when only one of many parameters is 
allowed to vary constitute decisive evidence against that theory. 

In our view, the issue of temperature dependence cannot be settled until experiments are able to detect a variation 
of the Casimir force with temperature. Clearly, such experiments are difficult. The precise experiments at very short 
distances (~ 100 nm), where the Casimir forces are largest, are not the best place to look for temperature variation, 
for the temperature dependence is relatively small there. Rather, experiments should be conducted at the micrometer 
scale, where the effects, if our theory is correct, are above the 10% level, and approach a factor of two reduction from 
the ideal metal limit for larger separations. (The Lamoreaux experiment [Tsj . conducted at the micrometer separation 
scale, was probably not as accurate as claimed.) Proposals to perform measurements of the force between a cylinder 
and a plane [2(| and between eccentric cylinders [2l| have advantages because the forces are stronger than between a 
sphere and a plane, yet the difficulties in maintaining parallelism is not so severe as with two plane surfaces. In this 
way the possiblitity of measuring the temperature dependence at relatively large distances where thermal effects are 
largest may be accomplished. 

VI. NEW CALCULATIONS 

To this end, we present some new calculations of the force between parallel plates based on our theory, compared 
with that of Ref. The numerical computations are done on the basis of Eq. (4.18) in our paper and follow 
the recent paper of Bentsen et al. |22(. The maximum value of the quantity y — \Jk\-\- Q n a, a being the plate 
separation, is chosen to be y max = 30. In all calculations we impose the tolerance for the integrals to be 10~ 12 . The 
overall tolerance in the calculated sum over Matsubara numbers is taken to be 10 -8 . This last tolerance determines 
the highest Matsubara frequencies occurring in the m sum. These values are shown in the tables. 

As for permittivities, we are using the (new) data received from Astrid Lambrccht. These data extend from about 
1.5 x 10 11 rad/s to about 1.5 x 10 18 rad/s. An important virtue of our tables is that they show which frequency 
domains we actually use. The largest frequency region naturally occurs when the temperature is low. Thus for T = 1 
K, the temperature that we associate with T = on physical grounds, we use the region from about 8 x 10 11 rad/s 
to about 2 x 10 16 rad/s. That means, we are, even at this low temperature, working with frequencies that lie entirely 
within the region of Lambrecht's data. This implies that we do not have to involve the Drude relation at all, for the 
finite frequencies. Of course, the Kramers-Kronig relation, required by causality, is employed. There is no analytical 
approximation involved in our formalism, for any finite frequency. The only exception is the zero frequency case. We 
then need the Drude relation to assure that there is no contribution to the force from the zero frequency TE mode. 
This zero frequency contribution to the Casimir force is found analytically, not numerically. 

For comparison, we show in Table Q] the results given in Ref. [3j for the Leontovich impedance and plasma models. 
These results are for zero temperature. In their approach the temperature dependence is negligible. It will be noted 
that their two models do not agree, even though the authors seem to imply that either model is equally good. Our view 
is that the difference between these two models may be taken as a rough gauge of the accuracy of their predictions. 
Our results, with details about the frequency range used and the number of terms in the Matsubara sum, are given 
in Tables ITT1 IIVI for T = 1 K (sufficiently close to T = K, but a temperature at which our numerical technique is 
stable), T = 300 K, and T = 350 K. The latter is given in the hope that the temperature variation over a 50 K range 
near room temperature may soon become accessible to experiment. 
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Separation Impedance Method Plasma Method 

nm Ref. [3] Ref. [3] 

160 1144 1114.9 

200 509.3 501.8 

250 224.7 223.1 

400 38.90 38.98 

500 16.70 16.76 

700 4.605 4.628 



TABLE I: Casimir pressure between parallel plates in the two models discussed in Ref. ||. Pressures are given in mPa. 



Separation Pressure Highest frequency Number of terms 



nm 


mPa 


rad/s 


in sum 


160 


1144 


2.112 x 10 io 


25674 


200 


508.2 


1.713 x 10 16 


20824 


250 


223.7 


1.388 x 10 16 


16869 


400 


38.61 


8.875 x 10 15 


10789 


500 


16.56 


7.168 x 10 15 


8714 


700 


4.556 


5.187 x 10 15 


6305 


1000 


1.143 


3.674 x 10 15 


4466 



TABLE II: Our results for the Casimir pressure between gold plates, when T — 1 K. The first nonvanishing Matsubara 
frequency (corresponding to m — 1) is 8.226 x 10 11 rad/s. The highest frequencies are dependent on the plate separation a, as 
shown, and are determined by our chosen tolerance 10 -8 for the sum. The last column gives the number of terms in the sum. 

VII. CONCLUSIONS 

Were the controversy over the temperature dependence a purely theoretical one, we believe that a consensus along 
the lines we have been advocating, justified by realistic model calculations by Jancovici and Samaj Q, Buenzli 
and Martin |15| . and Sernelius would now be nearly universally accepted. Evidently what has prevented this 
consensus from developing is the apparent conflict that the resulting predictions have with experiment. However, we 
believe it is premature to regard that conflict as decisive. It is imperative to perform room-temperature experiments at 
larger separations^since the relative temperature corrections are then much larger. The optimum separation seems to 
lie around 2 /im [23j . Indeed, the early experiments of Lamoreaux ^ij support the Drude theory at large separation, 
although the error bars are large. The experimental situation when a is about 1 /im or less (which is the main region 
of the Decca experiment may be influenced by extra factors coming into play at short distances, such as surface 
plasma effects 24]. At short distances, the results are moreover very sensitive to differences in separation; an error of 
1 nm may give rise to errors in the pressure of about 10 mPa or 2% at a = 200 nm. 
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Separation Pressure Highest frequency Number of terms 



nm 


mPa 


rad/s 


in sum 


160 


1127 


2.122 x 10 lb 


86 


200 


497.8 


1.727 x 10 16 


70 


250 


217.6 


1.407 x 10 16 


57 


400 


36.70 


8.884 x 10 15 


36 


500 


15.49 


7.403 x 10 15 


30 


700 


4.127 


5.429 x 10 15 


22 


1000 


0.9852 


3.702 x 10 15 


15 



TABLE III: Same as in Table El but at T = 300 K. The first nonvanishing Matsubara frequency is 2.468 x 10 14 rad/s. The 
highest frequencies are determined by the same tolerance 10 -8 for the m sum as before. 
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Separation Pressure Highest frequency Number of terms 



nm 


mPa 


rad/s 


in sum 


160 


1124 


2.131 x 10 lb 


74 


200 


495.7 


1.727 x 10 16 


60 


250 


216.4 


1.411 x 10 16 


49 


400 


36.35 


8.925 x 10 15 


31 


500 


15.30 


7.198 x 10 15 


25 


700 


4.052 


5.470 x 10 15 


19 


1000 


0.9590 


3.743 x 10 15 


13 



TABLE IV: Same as in Table EU but at T = 350 K. The first nonvanishing Matsubara frequency is 2.879 x 10 14 rad/s. The 
highest frequencies are determined by the same tolerance 10 -8 for the m sum as before. 
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